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ABSTRACT 

Type la supernovae (SNe la) play an important role in diverse areas of as- 
trophysics, from the chemical evolution of galaxies to observational cosmology. 
However, the nature of the progenitors of SNe la is still unclear. In this pa- 
per, according to a detailed binary population synthesis study, we obtained SN 
la birthrates and delay times from different progenitor models, and compared 
them with observations. We find that the Galactic SN la birthrate from the 
double-degenerate (DD) model is close to those inferred from observations, while 
the birthrate from the single-degenerate (SD) model accounts for only about 
1/2—2/3 of the observations. If a single starburst is assumed, the distribution of 
the delay times of SNe la from the SD model is a weak bimodality, where the WD 
+ He channel contributes to the SNe la with delay times shorter than 100 Myr, 
and the WD + MS and WD + RG channels to those with age longer than 1 Gyr. 

Subject headings: binaries: close — stars: evolution — supernovae: general — 
white dwarfs 



1. Introduction 

Type la supernovae (SNe la) appear to be good cosmological distance indicators due to 
their high luminosities and remarkable uniformity, and thus are used for determining cosmo- 
logical parameters(e.g., Riess et al. 1998; Perlmutter et al. 1999). They also throw light on 
the understanding of galactic chemical evolution owing to the main contribution of iron to 
their host galaxies. However, several key issues related to the nature of their progenitor sys- 
tems and the physics of the explosion mechanisms are still not well understood(Hillebrandt 
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& Niemeyer 2000; Wang et al. 2008a), and no SN la progenitor system before the explosion 
has been conclusively identified. This may raise doubts about the distance calibration, which 
is purely empirical and based on the SN la sample of the low red-shift universe. 

It is generally believed that SNe la are thermonuclear explosions of carbon-oxygen white 
dwarfs (CO WDs) in binaries. Over the past few decades, two famihes of SN la progenitor 
models have been proposed, i.e., the double-degenerate (DD) and single-degenerate (SD) 
models. The DD model involves the merger of two CO WDs, where the total mass of 
the two CO WDs is larger than the Chandrasekhar (Ch) mass limit (Webbink 1984; Han 
1998). For the SD model, the companion is probably a main-sequence (MS) star or a slightly 
evolved subgiant star (WD -|- MS channel), or a red giant star (WD -|- RG channel), or an 
He star (WD -|- He star channel) (Hachisu et al. 1999; Li & van den Heuvel 1997; Han & 
Podsiadlowski 2004; Chen & Li 2009; Meng et al. 2009; Lii et al. 2009; Wang et al. 2009a; 
Wang, Li & Han 2010). Note that, recent observations have suggested that at least some 
SNe la can be produced by a variety of progenitor systems (e.g., Patat et al. 2007; Wang et 
al. 2008b). 

At present, various progenitor models of SNe la can be examined by comparing the 
distribution of the delay time (between the star formation and SN la explosion) expected 
from a progenitor model with that of observations (e.g., Chen & Li 2007; Wang & Han 2010). 
There are three important observational results for SNe la, i.e., the strong enhancement of 
the SN la birthrate in radio-loud early-type galaxies, the strong dependence of the SN la 
birthrate on the colors of the host galaxies, and the evolution of the SN la birthrate with 
redshift. Mannucci et al. (2006) found that the observational results can be best matched 
by a bimodal delay time distribution, in which about half of the SNe la explode soon after 
starburst, with a delay time less than ~100 Myr, while those remaining have a much wider 
distribution with a delay time of ~3Gyr. It is suggested that 10% (weak bimodahty) to 
50% (strong bimodality) of all SNe la belong to the young SNe la (Mannucci et al. 2008). 
The existence of the young and old populations of SNe la is also supported by some other 
observations (e.g., Aubourg et al. 2008; Totani et al. 2008). 

The purpose of this paper is to study SN la birthrates and delay times from different 
progenitor models, and compare them with observations. We describe the binary population 
synthesis (BPS) approach for the different progenitor models in Section 2 and present the 
BPS results in Section 3. Section 4 ends the paper with a discussion. 
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2. Binary population synthesis 



In order to investigate SN la birthrates and delay times, we performed a series of Monte 
Carlo simulations in the BPS study. In each simulation, we followed the binary evolution 
via the rapid binary evolution code developed by Hurley et al. (2002). 



The progenitor of an SN la is a close WD binary, which has most likely emerged from 
the common envelope (CE) evolution of a giant binary. During the binary evolution, the 
primordial mass ratio (primary to secondary) is crucial for the mass transfer. If it is larger 
than a critical mass ratio, the mass transfer may be dynamically unstable and a CE 
develops. The mass ratio varies with the evolutionary state of the primary at the onset 
of the Roche lobe overflow (RLOF) (e.g., Han et al. 2002). In this study we adopt Qc — 4.0 
when the primary is in the main-sequence stage or Hertzsprung gap. This value is supported 
by detailed binary evolution studies (e.g., Han et al. 2000). If the primary is on the first 
giant branch or asymptotic giant branch stage, we use 



where Mf is the mass of the primary, is the core mass of the primary, and x = 

d In Ri /d In Mf is the mass-radius exponent of the primordial primary and varies with com- 
position. If the mass donor stars (primaries) are naked He giants, Qc — 0.748 based on eq. 
(1) (see Hurley et al. 2002 for details). 

When a CE forms, embedded in the CE is a 'new' binary consisting of the dense core 
of the primary and the secondary. Owing to frictional drag within the envelope, the orbit 
of the 'new' binary decays and a large part of the orbital energy released in the spiral-in 
process is injected into the envelope. The CE ejection is still an open problem. Here, we 
use the standard energy equations (Webbink 1984) to calculate the output of the CE phase. 
The CE is ejected if 



where A is a structure parameter that depends on the evolutionary stage of the mass donor 
star, Mdon is the mass of the donor. Mace is the mass of the accretor, a is the orbital 
separation, M^nv is the mass of the donor's envelope, i?don is the radius of the donor, and the 
indices i and f denote the initial and final values. The right side of the equation represents 



2.1. Common envelope in bineiry evolution 



?c = [1.67-a; + 2(-^)^]/2.13. 



(1) 
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the binding energy of the CE, while the left side shows the difference between the final and 
initial orbital energy, and cicc is the CE ejection efficiency (the fraction of the released orbital 
energy used to eject the CE). For this prescription of the CE ejection, there are two highly 
uncertain parameters (i.e., A and cicc)- As in previous studies, we combine a^^ and A into 
one free parameter a^eX (Wang et al. 2009b). 

2.2. The progenitor model of SNe la 

Theoretically, there is a consensus that an SN la is the explosion and complete disinte- 
gration of a CO WD that has a mass close to the Chandrasekhar limit mass Mch (1.378 Mq 
in this work). We adopt two popular models for the progenitors of SNe la: the DD and SD 
models. 



2.2.1. The DD model 

For the DD model, SNe la arise from the merging of two close CO WDs that have a 
combined mass larger than or equal to the Mch (Webbink 1984; Han 1998). Both WDs 
are brought together by gravitational wave (GW) radiation on a timescale tow (Landau & 
Lifshitz 1971), 

tGw(yr) = 8 X lO^(yr) x ^"^^ +^^^^''V ^/3(h), (3) 
where P is the orbital period in hours, ^gw in years and Mi, M2 in Mq. 

The delay time from the star formation to the occurrence of an SN la is equal to the 
sum of the timescale on which the secondary star becomes a WD and the orbital decay time 
tow- Here, we set the ctceA = 1.5, which reproduces the number of the double-degenerate 
objects in the Galaxy (Hurley et al. 2002). 

2.2.2. The SD model 

For the SD model, we considered the WD + He star, WD + MS and WD + RG channels 
in this paper. Take the WD -|- He star channel for example, we followed the evolution of 
1 X 10^ sample binaries from the star formation to the formation of the WD -|- He star 
systems. We assumed that, if the parameters of a CO WD -|- He star system at the onset 
of the RLOF are located in the SN la production regions in the (logP\ M2) plane (Figure 
8 of Wang et al. 2009a), where and are the orbital period and the mass of the He 
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companion star at the onset of the RLOF, respectively, an SN la is produced. Note that, 
the method of the BPS study for the WD + MS and WD + RG channels is similar to that 
of the WD + He star channel. For the SN la production regions of these two channels, see 
Wang, Li & Han (2010). Here, we set the ciccA = 0.5, which is our standard model for the 
formation of these SD channels (Wang, Li & Han 2010). 



2.3. Basic peirameters for Monte Ceirlo simulations 

In the BPS study, the Monte Carlo simulation requires as input the initial mass function 
(IMF) of the primary, the initial mass-ratio distribution, the distribution of initial orbital 
separations, the eccentricity distribution of binary orbit, and the star formation rate (SFR) 
(Wang et al. 2009b; Han 2008a,b). 

(1) The IMF of Miller & Scalo (1979) is adopted. The primordial primary is generated 
according to the formula of Eggleton et al. (1989), 

0.19X 

= (1 - X)0.75 + 0.032(1 - X)o.25^0' (4) 

where X is a random number uniformly distributed in the range [0, 1] and Mf is the mass 
of the primordial primary, ranging from 0.1 Mq to 100 Mq. 

(2) The initial mass-ratio distribution of the binaries, q', is quite uncertain for binary 
evolution. For simplicity, we take a constant mass-ratio distribution, 

n{q') = 1, < g' < 1, (5) 

where q' = M^/Mf. 

(3) We assume that all stars are members of binaries and that the distribution of sepa- 
rations is constant in log a for wide binaries, where a is separation and falls off smoothly at 
a small separation 

a . n{a) = { ^-Ma^^^^, « < «o, 

t Qisep, ao < a < ai, 

where ccgep ~ 0.07, Oq = 10 it!©, Oi = 5.75 x 10^ it!© = 0.13 pc and m 1.2. This distribution 
implies that the numbers of wide binaries per logarithmic interval are equal, and that about 
50 per cent of stellar systems have orbital periods less than 100 yr (Han et al. 1995). 

(4) A circular orbit is assumed for all binaries. The orbits of semidetached binaries 
are generally circularized by the tidal force on a timescale which is much smaller than the 
nuclear timescale. Also, a binary is expected to become circularized during the RLOF. 
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(5) We simply assume a constant SFR over the past 15 Gyr, or, alternatively, a delta 
function, i.e., a single starburst. In the case of the constant SFR, we calibrate the SFR by 
assuming that one binary with a primary more massive than 0.8 Mq is formed annually (see 
Han et al. 1995; Hurley et al. 2002). From this calibration, we can get SFR = 5MQyr~^. 
For the case of the single starburst, we assume a burst producing 10^^ Mq in stars. In fact, 
a galaxy may have a complicated star formation history. We choose only these two extremes 
for simplicity. A constant SFR is similar to the situation of spiral galaxies (e.g., Han & 
Podsiadlowski 2004), while a delta function to that of eUiptical galaxies or globular clusters. 



3. The results of bineiry population synthesis 

3.1. The birthrates of SNe la 

We performed four sets of simulations with metallicity Z — 0.02 and SFR = 5 MqYi'^ 
to systematically investigate Galactic birthrates of SNe la for the DD and SD models, where 
the SD model includes the WD + He star, WD + MS and WD + RG channels (see Figure 
1). We find the birthrate from the DD model ~2.9 x 10~^ yr~^, only slightly lower than 
the birthrate inferred from observations (i.e., 3—4 x lO^'^ yr~^; Cappellaro & Turatto 1997), 
while the total birthrates from the SD models can only account for about 1/2—2/3 of the 
observations, where the birthrate from the WD + He star channel ~0.3 x 10~^ yr~^, the 
WD + MS channel ~1.8 x 10"^ yr"^ and the WD + RG channel ~3 x 10"^ yr-^ 

The SN la birthrate in galaxies is the convolution of the distribution of the delay times 
(DDT) with the star formation history (SFH): 

v{t) = f SFR{t - t')DDT{t')dt', (7) 
Jo 

where the SFR is the star formation rate, and t' is the delay times of SNe la. Due to a 
constant SFR adopted in this paper, the SN la birthrate ^{t) is only related to the DDT, 
which can be expressed by 

r 0, t <ti, 

DDT{t) = <j DDT'{t), ti<t< ta, (8) 
I 0, t> t2, 

where ti and ta are the minimum and maximum delay times of SNe la, respectively, and the 
DDT' is the distribution of the delay times between ti and ta- If t is larger than ta, eq. (7) 
can be written as ^ 

u{t) = SFR / DDT\t')dt' = constant. (9) 
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Therefore, the SN la birthrates shown in Figure 1 seem to be completely flat after the first 
rise. 



3.2. The delay times of SNe la 

Figure 2 displays the evolution of SN la birthrates for a single starburst with a total 
mass of 10^^ Mq. In the figure we see that SNe la from the DD model have the delay times 
of ~89Myr— 15Gyr, which are close to the observational results from Totani et al. (2008). 
However, it is suggested that the DD model is likely to lead to an accretion-induced collapse 
rather than to an SN la (Nomoto & Iben 1985). Thus, the DD model is not supported 
theoretically. We also find that SNe la from the SD models have a wide distribution of 
the delay times, where the WD + He star channel contributes to the SNe la with delay 
times shorter than 100 Myr, and the WD + MS and WD + RG channels to those with age 
longer than 1 Gyr (the WD + MS channel also contributes to the SNe la with intermediate 
delay times ~100Myr-l Gyr). Note that, Chen & Li (2007) studied the WD + MS channel 
by considering a circumbinary disk which extracts the orbital angular momentum from the 
binary through tidal torques. This study also provides a possible way to produce such old 
SNe la (~l-3Gyr). 

The SD model is currently a favorable progenitor model of SNe la. The distribution of 
the delay time from the SD models is similar to that derived from observations by Mannucci et 
al. (2006), except that the peak value of young population is smaller than that in Mannucci 
et al. (2006). The WD + He star channel produces 14 per cent of all SNe la, which 
constitutes the weak bimodality as suggested by Mannucci et al. (2008). The results from 
the SD models seem to be slightly smaller than those in Totani et al. (2008). However, given 
the high possibility of errors in the observation by Totani et al. (2008), our results are close 
to observations. 



4. Discussion 

In our BPS studies, we assumed that all stars are in binaries and about 50 per cent 
of stellar systems have orbital periods less than 100 yr. In fact, this is known to be a 
simplification. The binary fractions may depend on metallicity, environment and spectral 
type. If we adopt 40 per cent of stellar systems with orbital periods below 100 yr by adjusting 
the parameters in eq. (6), we estimate that the Galactic SN la birthrate from the DD model 
will decrease to be ~2.3 x 10"^ yr-^ the WD + He star channel ~0.24 x 10"^ yr-^ the WD 
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+ MS channel ~1.4 x 10"^ yr"^ and the WD + RG channel ~2.4 x 10"^ yp-^ 

Hachisu et al. (2008) investigated new evolutionary models for SN la progenitors, in- 
troducing the mass-stripping effect on a MS or a slightly evolved companion star by winds 
from a mass-accreting WD. Though the model offers a possible way to produce young SNe 
la, the model depends on the efficiency of the mass-stripping effect. We also find that the 
model produces very few young SNe la by a detailed BPS approach. Thus, we consider the 
WD + He star channel as a main contribution to the young population of SNe la, possible 
the whole young population. 

The Galactic SN la birthrate from the WD + RG channel is ~3xl0~^ yr~^, which is low 
compared with observations. Thus, SNe la from this channel may be rare. However, further 
studies on this channel are necessary, since this channel may explain some SNe la with long 
delay times. RS Oph and T CrB, both recurrent novae are probable SN la progenitors and 
belong to the WD -|- RG channel (e.g. Belczynski & Mikolajewska 1998; Hachisu et al. 
1999, 2007). Detecting Na I absorption lines with low expansion velocities indicates that 
the companion of the progenitor of SN 2006X may be an early RG star (Patat et al. 2007). 
We note that a symbiotic star with aspherical stellar wind may also provide a way for the 
evolution of the WD -|- RG systems towards SNe la (Lii et al. 2009). 

Provided that the DD model can produce SNe la, an explosion following the merger of 
two WDs would leave no remnant, while the companion star in the SD model would survive 
and potentially be identifiable (Wang & Han 2009). Thus, it will be a promising method to 
test different SD models of SNe la by identifying their surviving companions. 

The young and old populations of SNe la may have an effect on models of galactic 
chemical evolution, since they would return large amounts of iron to the interstellar medium 
much earlier or later than previously thought. In future investigations, we will explore the 
detailed influence of SNe la on the chemical evolution of stellar populations. 

This work is supported by the National Natural Science Foundation of Ghina (Grant 
No. 10821061), the National Basic Research Program of China (Grant No. 2007CB815406), 
and the Yunnan Natural Science Foundation (Grant No. 08YJ041001). 
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Fig. 1. — Evolution of Galactic SN la birthrates for a constant star formation rate {Z — 0.02, 
SFR = 5 Mgyr"^). The key to the line-styles representing different progenitor models is given 
in the upper left corner. 
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Fig. 2. — Similar to Figure 1, but for a single starburst with a total mass of lO^^M©. The 
observational points are from Totani et al. (2008) '^^1 



